Driven vortices in 3D layered superconductors: Dynamical ordering along the c-axis 



Alejandro B. Kolton/ Daniel Domfnguez,^ Cynthia J. Olson^ and Niels Gr0nbech- Jensen'^'''' 
^ Centra Atomico Bariloche, 8400 S. C. de Bariloche, Rio Negro, Argentina 
^ Department of Physics, University of California, Davis, California 95616 
^ Department of Applied Science, University of California, Davis, California 95616 
* NERSC, Lawrence Berkeley National Laboratory, Berkeley, California 94720 

(February 6, 2008) 

We study a 3D model of driven vortices in weakly coupled layered superconductors with strong 
pinning. Above the critical force Fc, we find a plastic fiow regime in which pancakes in different layers 
are uncoupled, corresponding to a pancake gas. At a higher F, there is an "smectic fiow" regime 
with short-range interlayer order, corresponding to an entangled line liquid. Later, the transverse 
displacements freeze and vortices become correlated along the c-axis, resulting in a transverse solid. 
Finally, at a force Fs the longitudinal displacements freeze and we find a coherent solid oi rigid lines. 

PACS numbers; 74.60. Ge, 74.40. +k, 05.70.Fh 



It is well-known that an external current can induce 
an ordering of the vortex structure in superconductors 
with pinning For a long time, it was believed that 
the high-current phase would have crystalline order. Re- 
cently, it has been found that different kinds of order are 
possible at high currents, depending on pinning strength 
and dimensionality |^-|^ . This has led to numerous theo- 
retical [^U^ , experimental and numerical studies . 
A crystal-like structure, which could be either a perfect 
crystal [|j or a Bragg glass is only possible in d = 3 
at large drives. In d = 2, or in d = 3 for intermedi- 
ate currents, a transverse glass is expected, with order 
only in the direction perpendicular to the driving force 
In the equilibrium vortex phase diagram, the 
behavior of vortex line correlations along the direction 
of the magnetic field (c-axis) has been intensively dis- 
cussed both experimentally Q and theoretically In 
the case of driven vortices, little is known on how the 
c-axis line correlations would behave in the different dy- 
namical regimes. Here we will address this issue starting 
from the less favorable case: weakly coupled supercon- 
ducting planes with strong pinning. We will show how 
the order along the c-axis and the in-plane structural or- 
der take place in a sequence of dynamical phases upon 
increasing current. 

We study pancake vortices in a layered superconduc- 
tor, considering the long-range magnetic interactions be- 
tween all the pancakes and neglecting Josephson coupling 
This model is adequate when the interlayer peri- 
odicity d is much smaller than the in-plane penetration 
length A|| jl0| . Previous simulations of driven vortices in 
3D superconductors have been performed using Langevin 
dynamics of short-range interacting particles [Tl| or the 
driven isotropic 3D XY model [ p^ . 

The equation of motion for a pancake located in posi- 
tion Ri = {Vi, Zi) = {xi,yi,nid), (z ee c), is: 



(1) 



where pij = |ri — Vj\ and — \zi — Zj\ are the in- 
plane and inter-plane distance between pancakes 
Pip = |ri — Tpl is the in-plane distance between the vortex 
i and a pinning site at Rp — (rp,Zi), ry is the Bardeen- 
Stephen friction, and F = -^J x z is the driving force 
due to an in- plane current J. We consider a random 
uniform distribution of attractive pinning centers in each 
layer with Fp = —2Ape^^P^°-p'' r/a^, where Op is the pin- 
ning range. The magnetic interaction between pancakes 
F^p, z) = Fp{p, z)f is given by 0,||: 



P 

Fp(p,Zn) = -J^ — 



I-^fl-e-^A" 
A V 



(2) 
(3) 



Here, R = ^/z^Tp^ and A = 2A|j/d is the 2D thin- 
film screening length. An analogous model to Eqs. (2- 
3) was used in |l4|. We normahze length scales by A||, 
energy scales by Ay — (fy^/An'^A, and time is normal- 
ized by T = rjX^^/Ay. We consider iV^, pancake vortices 
and Np pinning centers per layer in TV; rectangular lay- 
ers of size Lx X Ly, and the normalized vortex density is 
Uy ^ B Xy^o ^ (ao/A||)^. We consider Uy = 0.29 with 

Ly = 16A|| and = V^/2Ly, Ni ^ 8 and Ny 64. 
We take a pinning range of ap — 0.2, a large pinning 
strengh of Ap/Ay =0.2, with a high density of pinning 
centers Up = 3.125ni,. The model of Eq.(2-3) is valid in 



the hmit d < A|| < A. We take d/Xu 
corresponds to BSCCO compounds llQl . 



= O.OI, which 
Moving pan- 
cake vortices induce a total electric field E = — v x z, 
with V = N Ni Si study the dynamical regimes 

in the velocity-force curve at T = 0, solving Eq. (1) for 
increasing values of F = Fy . We use periodic bound- 
ary conditions both in the planes and in the z direction 
and interactions between all pancakes in all layers are 
considered ]T^ . The periodic long-range in-plane and 
inter-plane interaction is evaluated using Ref. JlSf . The 
equations are integrated with a time step of At = O.OIr 
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and averages are evaluated in 16384 integration steps af- 
ter 2000 iterations for equilibration. Each simulation is 




FIG. 1. Vortex trajectories in the first five layers: (a) 
F=0.6, (b) F=l.l, (c) F=2.0. (d) F=3.9. Surface intensity 
plot of the averaged in-plane structure factor S(k): (e) F=0.6, 
(f) F=l.l, (g) F=2.0. (h) F=3.9. 

started at F = with a triangular vortex lattice and 
slowly increasing the force in steps of AF — 0.1 up to 
values as high as F = 8. 

We start with a qualitative description of the differ- 
ent steady states that arise as a function of increas- 
ing force. In Figure l(a-d) we show the vortex trajec- 
tories {Ri(<:)} for typical values of F by plotting the 
positions of the pancakes in five of the layers for all 
t. In Fig.l(e-h) we show the average in-plane struc- 
ture factor S{k) = |^ exp[ik.r„(t)]|2), with 

k = {kx,ky). Above the depinning critical force Fc, 
we find the following dynamical phases. (i)Plastic flow 
{Fc < F < Fp): Pancakes flow in an intrincate network 
of "plastic" channels similar to the behavior found in 
2D 1^,0. The motion in different planes is completely 
uncorrelated, [Fig. 1(a)] and there is no signature of or- 
der in the structure factor [Fig. 1(e)]. (ii) Smectic flow 
{Fp < F < Ft): The motion organizes in "elastic" chan- 
nels that are almost parallel and separated by a distance 
~ ap, see Fig. 1(b). Small and broad "smectic" peaks 
appear in ^(k) for k • F = [Fig. 1(f)]. There are "acti- 
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FIG. 2. (a) Velocity-force curve, left scale, black points, 
dV/dF (differential resistance), right scale, white points, (b) 
Intensity of the Bragg peaks. For smectic ordering S{Gi), 
Ky = 0, (x) symbols. For longitudinal ordering S{G2,3), 
7^ 0, {+) symbols, (c) Diffusion coefficient for trans- 
verse motion D^c, (A), left scale. Longitudinal displacements 
{[Ay{t)]'^) for a given t as a function of F, (■), right scale. 

each other between neighboring planes, (iii) Transverse 
solid {Ft < F < Fs): There are well defined channels 
in all the planes and the pancakes do not jump between 
channels [Fig. 1(c)]. The structure factor has sharp smec- 
tic peaks and small "longitudinal" peaks (k • F ^ 0) 
have appeared [Fig. 1(g)]. The location of channels is 
correlated in the c-axis. {iy)Coherent solid {F > Fg): 
The channels become more straigth with small transverse 
wandering [Fig. 1(d)]. The S'(k) shows well defined peaks 
for all k in the reciprocal lattice [Fig. 1(h)]. The dynami- 
cal phases (i)-(iii) are similar to the ones found previously 
in 2D thin films §. 

The characteristic forces Fc, Fp, Ft, Fg separating the 
different dynamical phases are obtained from the analysis 
of the in-plane and out of plane structural and dynamical 
correlations. We show in Fig. 2 the in-plane structure fac- 
tor and temporal fluctuations, obtained in the same way 
as for 2D In Fig. 2(a) we plot the average velocity 
V = {Vy{t)) = Y.n. in the direction of the 



force as a function of F and its corresponding derivative 
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dV/dF (differential resistance) . The force Fp corresponds 
to the peak in the differential resistance. We also see a 
small second maximum in dV/dF for a force between 
Ft and Fs In Fig. 2(b) we plot the magnitude of 

the peaks in the in-plane structure factor. We show the 
peak height at Gi = 27r/aox, corresponding to smectic 
ordering, and the average of the peaks corresponding to 
longitudinal ordering at G2 = ±27r/ao(l/2, V3/2) and 
G3 = ±27r/ao(-l/2, \/3/2). We see that at Fp the smec- 
tic peak rises up from zero, then at Ft it reaches an al- 
most constant value and later at Fg it has a small jump. 
The longitudinal peak has a small finite value for forces 
above Fp, and only at Fg shows a significant increment 
towards a large value. Comparing with the previous 2D 
results 0, we can make the reasonable assumption that 
for Fp < F < Ft there is only short-range smectic or- 
der (since there is activated transverse diffusion beteween 
elastic channels, see below), for Ft < F < Fg there is 
probably quasi-long range smectic order but short range 
longitudinal order, and above Fg there is both transversal 
and longitudinal order (quasilong-range or long-range). 
What is new, compared with the 2D thin film case 0, 
is that above a force Fg there is a significant amount of 
longitudinal order. This may correspond either to a mov- 
ing crystalline phase (if there is long-range order) or to 
a moving Bragg glass (if there is quasilong-range order) 
1^. We have verified that, for a given F > Fg, there is 
both longitudinal and transversal order for system sizes 
of AT/ X = 5 X 36, 5 X 64, 8 x 64, 8 x 100, 10 x 100. How- 
ever, a detailed finite size analysis is not possible with 
these few small samples. We complement our discussion 
of the in-plane physics with the study of the temporal 
fluctuations, which are shown in Fig. 2(c). We calcu- 
late the transverse diffusion coefficient from the aver- 
age quadratic transverse displacements of vortices from 
their center of mass position (X„,y"„), ^ 

XnM - a^»(0) + ^nM? ~ DJ. We find that is 
maximum at Fp in coincidence with the peak in dV/dF. 
Below Fp diffusion is through the intrincate network of 
plastic channels, above Fp diffusion is through activated 
jumps between elastic channels. sharply drops to 
zero at Ft, indicating that transverse displacements are 
localized in the transverse solid phase |^]. The drift 
from the center of mass of longitudinal displacements 
([A2/(t)]2) = ([y,(t) - y„,(t) - y,;(0) + Y„M?) is su- 
perdiffusivc for F < Fg, similar to the results observed in 
2D films 0. For F > Fg the longitudinal displacements 
become frozen in a constant value {[Ay{t)]'^) < gq/Ni, as 
it is shown in Fig. 2(c). Since in-plane displacements are 
localized and there are large transversal and longitudinal 
Bragg peaks, we call this phase a coherent solid. 

Let us now discuss how the ordering along the c-axis 
takes place. We analyze the pair distribution function: 
giP,n) = ^^{T.i^jS{p ~ Ptj)5n,n,,)- From 5r(p,n) 
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we define a correlation function along c-axis Cziji) 



FIG. 3. (a) Correlation parameter in c-direction of instan- 
taneous configurations C'z{n) vs F for n — 1, 2, 3, 4 interplane 
distance, (b) Trajectories overlap correlation parameter in 
c-direction On vs F for n = 1, 2, 3, 4 interplane distance, (c) 
Voltage fluctuations in c-direction (S^Vc) vs F. 

finite Cz{n = 1), meaning that pancakes in neighbor- 
ing planes are coupled and a "vortex line" can therefore 
be defined. In principle, an exponential decay Cz(n) 
exp(— n/fz) would define a correlation length for the vor- 
tex line . On the other hand, long-range ordering will 
be given by Cz(n 00) > 0. In Fig. 3(a) we 

show Cz{n) as a function of F for n ~ 1,2, 3, 4. We see 
that at Fp there is an onset of short-range order along 
the c-axis with a finite Cz{n — 1). At higher forces be- 
tween Fp and Ft the other Cz{n > 1) start to rise. The 
absence of correlations for F < Fp means that pancake 
motion is completely random and uncorrelated between 
different planes. Therefore, we propose that the plastic 
ffow regime corresponds to a pancake gas. Above Fp, in 
the smectic flow regime, it is possible to define a vortex 
line with short range correlations along the c-axis. Since 
there are in-plane jumps between elastic channels (i.e., 
cutting and reconnection of flux lines) we may consider 
this phase as an entangled line liquid. Above Ft, Cz{n) 
is finite for all n considered and tends to saturate upon 
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increasing n. This indicates that vortex Unas become 
more stiff above Ft. We also analyzed the c-axis corre- 
lation between averaged vortex densities. We first define 
/9„(r, n,t) = -w^J^i^i^- ^ni(t)) taking a coarse-graining 
scale Ar — ao/2 (results do not vary much for Ar — 
ao/4). The regions where the average density (yOu(r, n)) is 
large define the paths of steady state vortex motion. We 
can thereby calculate the overlap function of vortex tra- 
jectories between different planes as 0„ — Cp{n)/Cp{0), 
withCp(n) = h^[J2^ J dr{p,{r,m)){p,{r,m+n))]-l. 
This is shown in Fig. 3(b). We see that 0„ also has an 
onset at Fp. For Fp < F < Ft, we have some overlap 
of the elastic channels that decreases with increasing n, 
consistent with the entangled line-liquid picture. More 
interestingly, at Ft the overlap function 0„ becomes inde- 
pendent of n. This means that there is long-range c-axis 
coupling of the path of the elastic channels. When trans- 
verse displacements become localized in the x-direction, 
they also become locked in the c-direction. Thus, the 
freezing of in-plane transverse displacements occurs si- 
multaneously with a transverse disentanglement of flux 
lines at Ft. A striking result is that we find 0„ ~ 1 
above F^, i.e., a perfect c-axis coupling of elastic channels 
(within the scale ~ ao/4). Another interesting point to 
consider is the correlation of vortex velocities. If vortices 
in different planes move at different velocities, they will 
induce a Joscphson voltage difference along the c-axis 
given by Vn,n+i{r,t) = ^^0„,„+i(r, i), with 0„,„+i 
the superconducting phase difference between planes n 
and n -|- 1. A good approximation for pancakes at r„^i is 
to write (j)n,n+i{r,t) = Y.i[fi'^-'^n,i)~f{r~r,i+i,i)] with 
/(r) aTctan{x/y). We can therefore estimate the c-axis 
voltage fluctuations as {S^Vc) = J2n I ii^n n+ii'^^t)) - 
(K,„+i (r, mdr^A E„ [{VI) - {Vnf] ~ [(V„ • V„+i) - 
(V„) • (V„+i)]; with V„(i) = ^ J2i Vn,i(i), and the con- 
stant A ^ log A if L > A or A ^ log(i) otherwise. It 
is clear that (S'^Vc) = for pancakes moving with the 
same velocity in all planes. We see in Fig. 3(c) that the 
voltage fluctuations have a maximum at Fp. For F > Fp, 
{6'^Vc) decreases, and above Fg it reaches an almost F- 
independent value. The fact that (S'^Vc) does not vanish 
above Fg is consistent with the result that Cz{n) < 1 for 
all values of F in Fig. 3(a). In other words, while trans- 
verse displacements are strongly correlated along the c 
direction for large forces [Fig. 3(b)], the longitudinal dis- 
placements in different planes are weakly correlated . 

In conclusion, we have clearly distinguished different 
dynamical phases in 3D layered superconductors consid- 
ering both in-plane and c-axis ordering The onset 
of short-range c-axis correlations could be studied ex- 
perimentally with plasma resonance measurements [ pT| . 
The long-range ordering along the c-axis could be stud- 
ied through simultaneous measurements of pc resistivity 
and in-plane current- voltage response . 
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